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I . INTRODUCTION 


Observation  of  so-called  "precursor"  shocks  occurred  early- 
more  than  two  decades  ago-  in  nuclear  testing  programs.  In-the-field 
measurements  for  air  bursts  over  terrain  produced  characteristic 
pressure  traces  which  portray  (see  Figure  1)  a significant  pressure 
rise  leading  the  main  pressure  wave  arrival.  1 The  importance  of  the 
precursor  phenomenon  for  nuclear  weapons  lies  in  its  alteration  of  the 
blast  intensity  and  wave  form:  the  dynamic  pressure  is  signifi- 
cantly increased,  thereby  increasing  the  chances  for  defeat  of  struc- 
tures strengthened  against  static  pressure  forces  but  still  vulnerable 
to  drag  forces. 

The  mechanism  inducing  the  precursor  formation  in  the  field  is 
generally  believed  to  be  the  presence  prior  to  shock  arrival  of  a 
thermal  layer  adjacent  to  the  earth's  surface.  Such  a layer  may  be 
formed  by  irradiation  of  the  earth's  surface  by  the  nuclear  burst  and 
subsequent  evaporation  of  the  surface  into  the  adjoining  air  layer. 
This  heated-air  layer,  the  "thermal  layer"  so-called,  carries  with  it 
a higher  sound  speed,  with  the  important  result  that  waves  incident  on 
and  in  the  layer  may  travel  faster  in  it  and  thus  lead  the  portion  of 
incident  wave  traveling  in  the  colder  region,  thereby  giving  rise  to 
the  precursor. 


Field  measurements  of  the  air  temperature  prior  to  shock  arrival 
from  the  nuclear  burst  indicate  rather  high  temperatures,  T > 1200°K 

having  been  recorded  within  a few  feet  of  ground  at  ground  ranges 
GR  ^ 300-  600m. ^ The  temperature  history  as  well  as  its  variation 
with  height  is  less  clearly  known  from  these  experiments,  due  to 
difficulties  with  instrumentation  in  the  extreme  environment  and  with 
interpretation  of  records. 


1.  K.  Kaplan , P . J.  Norris,  P.  A.  Ellis,  and  D.  C.  Sachs,  " Nuclear 
Blast  Phenomena  - Vol.  II."  DASA  1200  (1971).  CONFIDENTIAL.  This 
reference  contains  an  extensive  bibliography  and  discussion  of 
theoretical  and  experimental  work  on  the  precursor. 

2.  T.  Ft.  Broida , A.  Broido,  and  A.  B.  Willoughby,  "Air  Temperature 
in  the  Vicinity  of  c Nuclear  Detonation"  Operation  TUMBLER 
Project  8.  2WT542. 

3..  E.C.Y.  Inn,  "Air  Temperature  Measurements  Over  Several  Surfaces." 
Operation  TEAPOT  Project  8.4e  WT-1149  (1957). 
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In  the  laboratory  at  about  same  time  period  the  precursor 
phenomenon  was  studied  by  Griffith**  and  by  Varwig  and  Zemel^  using 
plane  or  reflected  shocks  interacting  with  air  over  a heated  plate. 

Interesting  features  of  the  interaction  are  given,  and  Griffith  makes  com- 
parisons with  his  theoretical  model  derived  for  moderate  strength  shocks 
and  plate  temperatures.  Jahn^1  investigated  the  interaction  of  a plane 
shock  incident  on  two  gases  of  different  sound  speed  separated  by  an 
interface.  Discussion  of  significant  features  is  based  on  detailed 
knowledge  of  gasdynamics  and  is  invaluable  in  this  respect. 


Early  theoretical  treatments  of  a small  disturbance  incident 
on  the  interface  between  super-  and  subsonic  regions  and  subsequent 
propagation  of  the  disturbance  have  been  made  by  Howarth^  and  by  Tsien 
and  Finston®,  whose  interests  then  were  pointed  toward  shock-boundary 
layer  interactions.  The  principal  finding  appears  to  be  that  the 
disturbance  in  the  subsonic  region  propagates  upstream  and  ahead  of 
that  in  the  supersonic  region.  No  detailed  treatments  exist  insofar 
as  is  known  for  shocks  or  large  disturbances  analogous  to  the  linear- 
ized treatments  mentioned  above.  Simple  phenomenological  models*  of 
idealized  shock-thermal  layer  interaction  have  been  constructed  and 
analyzed  in  accordance  with  acoustical  theory  as  well  as  with  oblique 
shock  theory.  A major  prediction,  verified  by  laboratory  experiments, 
is  a relation  between  the  precursor  shock  angle  and  temperature  ratio 
of  (uniform)  hot  and  cold  regions: 


sin  p = (Tc/Th)1/2 


Also,  numerical  techniques  and  models  for  predicting  flows  of  great 


4 . W.  C.  Griffith 3 "Interaction  of  a Shock  Wave  with  a Thermal  Boundary 
Layer."  J.  Aero.  Sci.  23,  16-23  (1956). 

5.  R.  Varwig  and  J.  Z&nelj  "The  Interaction  of  Shock  Waves  with  a 
Thermal  Layer"  NAVORD  Report  4021 , AFSWR  266  (1955). 

6.  R.  G.  Jahn3  "The  Refraction  of  Shock  Waves  at  a Gaseous  Interface  - 
III.  Irregular  Refraction."  Princeton  University  Tech.  Report  1119 
(1955). 

7.  L.  Howarthj  "The  Propagation  of  Steady  Disturbances  in  a Supersonic 
Stream  Bounded  on  One  Side  by  a Parallel  Subsonic  Stream. " Proc. 
Combr.  Phil.  Soc.  44,  380-90  (1948). 

8.  H.  S.  Tsien  and  M.  Finston3  ". Interaction  Between  Parallel  Streams 
of  Subsonic  and  Supersonic  Velocities."  J.  Aero.  Sci.  163  515-528 
(1949). 
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complexity  have  been  and  are  being  developed. 

The  purpose  of  this  work  was  to  produce  a clean-gas,  model  thermal 
layer  with  precisely  known  conditions  and  to  study  the  interaction  of 
a plane  shock  with  the  thermal  layer.  Observation  by  means  of  a Mach- 
Zehnder  interferometer  would  allow  deducing  temperature  profiles 
existing  in  the  pre-shock  and  post -shock  regions  of  the  flow.  Results 
would  be  useful  for  guidance  and  comparison  with  numerical  models  and 
predictions.  The  layer  was  created  by  a heated  plate  set  in  the  test 
section  of  a shock  tube,  the  same  technique  as  used  by  Griffith4  at 
much  lower  temperatures.  A single  shock  condition  was  most  thoroughly 
investigated,  for  a range  of  plate  temperatures,  although  results  at 
one  or  two  other  conditions  were  obtained.  Results  for  shocks  of  small 
but  realistic  overpressures  (2.8-21)104  Pa  (4-30  psi)  were  thought 
desirable,  -with  thermal  layer  temperatures  of  < 1000°K.  An  accompa- 
nying objective  during  this  project  was  to  establish  feasibility  for 
methods  of  obtaining  a thermal  layer  which  would  allow  scaling  up  to 
a much  larger  area  thermal  layer  over  which  a small  explosive  could  be 
detonated,  the  rationale  here  being  that  the  computer  codes  existing 
could  model  explosives  and  hence  a more  severe  test  of  codes  would  be 
afforded. 


II.  EXPERIMENTAL  APPARATUS  AND  PROCEDURES 
A.  Shock  Generation  and  Conditions 


As  mentioned  above,  a plane  shock  generated  in  a shock  tube 
interacts  with  a thermal  layer  created  by  a heated  plate  mounted  in 
the  test  section.  An  over-all  sketch  of  the  experimental  setup  is 
shown  in  the  Figure  2.  The  desired  shock  overpressures  were  obtained 
approximately  by  use  of  driver-to-driven  pressure  ratios  given  by  one- 
dimensional shock  tube  theory.  The  precise  flow  condition  at  the  test 
section  is  then  obtained  from  the  monitored  shock  velocity  by  use  of 
shock  tube  tables.  Fast  response,  commercially  available  quartz 
piezoelectric  pressure  gages  monitored  the  shock  velocity;  and  mea- 
surements between  gage  oscilloscope  traces  are  expected  to  give  shock 
velocity  to  within  ± 1%,  as  is  usual  in  shock  tube  work. 

B.  Mach-Zehnder  Interferometer 


The  Mach-Zehnder  interferometer  used  in  these  experiments  has 
25  cm  diameter  plates  recently  resurfaced,  to  adequately  image  the 
fine  fringes  seen  in  the  thermal  layer.  Backlighting  for  the  inter- 
ferometer was  provided  by  a 2J  passive-dye  Q-switched  ruby  laser. 

Aside  from  apparent  operation  under  Murphy's  Laws*,  the  laser  back- 

*Murphy’s  propositions  are,  according  to  a f.  n.,  p.  393,  H.  T.  Davis , 
Introduction  to  Nonlinear  Differential  and  Integral  Equations:  1) 

If  anything  can  go  wrong , it  will ; 2)  Things  when  left  alone  can  only 
go  from  had  to  worse ; 3)  Nature  sides  with  the  hidden  flaw. 
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lighting  provided  ample  intensity  and  a well  defined  wavelength  for 
the  data  reduction  procedures. 


a problem  with  graininess  in  the  film  image  was  noted  and  caused 
some  concern.  A disc  of  ordinary  sawn  quartz  had  been  in  use^  to 
diffuse  and  scatter  adequately,  thus  to  broaden  and  smooth  out  intensity 
variations  in  the  small  laser  (<  1 mm  diam)  beam  before  its  entering  the 
collimating  lens.  However,  repeated  firings  at  the  high  power  densities 
had  apparently  caused  deterioration  of  the  special  properties  inherent 
in  the  particular  area  used  of  the  quartz  disc,  resulting  in  objec- 
tionable fixm  graininess.  And  efforts  to  find  another  suitable  spot 
on  this  and  other  such  discs  were  abandoned,  when  chances  for  success 
appeared  minimal  and  not  clearly  repeatable. 


An  obvious  beam  broadener  is  the  ordinary  ground  glass  screen. 
However,  observation  in  uhe  film  plane  with  a small  cw  laser  shows 
merely  that  granularity  associated  with  ground  glass.  But,  motion 
of  the  glass  in  its  plane  (normal  to  optical  axis)  produces  an 
acceptably  evenly  illuminated  field,  an  observation  noted  several 
years  earlier. However,  such  motion  and  transport  during  pulsed 
operation  are  beyond  present  technology,  and  this  sort  of  modification 
to  the  optical  train  was  not  entertained.  Instead,  it  was  hoped  that 
an  effect  analogous  to  plate  motion  could  be  achieved  by  simply  u:ing 
two  ground  glass  plates  in  the  beam,  with  random  granular  structuies 
of  the  two  plates  falling  again  randomly  in  the  film  plane,  but  now 
more  evenly  distributed.  Happily  such  was  the  case,  and  the  film 
plane  was  once  again  satisfactorily  illuminated.  Trials  showed  that 
other  ground  glass  pairs  could  be  thrust  into  the  beam  with  no  apparent 
loss  in  quality  of  illumination.  A small,  uncritical  separation  between 
plates  was  seemingly  necessary  to  obtain  optimum  contrast.  All  of  the 
pictures  produced  in  this  report  were  made  with  this  simple  beam 
broadener. 


C.  Heated  Plate  Thermal  Layer 

An  important  goal  in  this  project  was  to  devise  a means  of  pro- 
duction of  a regular,  well-defined  thermal  layer.  The  technique 
chosen  would  lend  itself  as  mentioned  in  the  Introduction,  to  scaling 
upward  for  producing  a similar  layer  of  much  larger  extent  for  testing 
with  explosives.  A number  of  methods  was  considered:  1)  flowing  of 

a heated  layer  of  air  over  test  plate  or  test  section  wall;  2)  rapid 
heating  of  the  air  above  a plate  via  a grid  of  suddenly  energized 


9.  F.  H.  Oertely  "Laser  Interferometry  of  Unsteady  Underexpanded  Jets." 
BEL  E 1694  (1974)  AD  773664. 

10.  T.  A.  Wiggins 3 "Methods  for  Spreading  Laser  Beams.  " Appl.  Optics  10, 
963-5  (1971). 
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electrical  wires  or  ribbon;  3)  electrical  resistance  heating  of 
a plate  which  then  acts  to  heat  the  adjacent  air  at  the  testing 
station;  4)  heating  of  a suitable  plate  via  microwaves,  possibly 
making  use  of  surplus  radar  antennas  (?);  5)  burning  of  a layer  of 
gaseous  fuel.  Because  of  experimental  limitations  the  decision  was 
made  to  proceed  with  3) , the  brute-force  method,  on  the  thought  that 
success  was  more  assured  both  in  producing  the  thermal  layer  at  desired 
temperatures  and  in  scaling  upward  to  the  larger  thermal  layers. 

Results  of  the  homemade  electrical  furnace  are  shown  in  the  Figure  3. 
Only  construction-type  fire  brick  (thermal  conductivity  ^0.6  Btu/hr. 
ft.°F)  was  used  since  on  hand,  aid  these  were  laid  into  top  and 
bottom  sections  of  the  piece  of  30  cm  diameter  aluminum  pipe  to  form 
a heating  cavity  for  the  41  X 10  X 1.25  cm  steel  test  plate.  Resistance 
wire  was  used  to  allow  'v  1000  W heating  rate  as  an  initial  trial. 

The  commercial  availability  of  much  better  insulating  materials  in 
both  brick  and  wool  blanket  form  (thermal  conductivity  0.03  Btu/hr. 
ft.  F)  suggest  much  greater  heating  efficiencies  are  possible  over  the 
present  crude  assembly.  Thus,  scaling  the  furnace  to  larger  sizes 
entails  simply  stringing  of  more  lengths  of  resistance  wire  and 
blanketing  with  more  insulating  material. 

A heating  rate  curve  for  the  oven,  as  measured  by  a thermocouple 
located  near  a wall  of  the  cavity  with  test  plate  in  place,  is  shown 
in  Figure  4a.  The  accompanying  heating  rate  for  the  test  plate,  as 
measured  by  a thermocouple  brazed  to  test  plate's  bottom  surface, is 
also  shown.  A cooling  curve  for  the  test  plate  mounted  in  the  test 
section  is  shown  in  Figure  4b.  The  thermocouple  readings  are  expected 
to  be  merely  indicative  of  the  plate's  test-surface  temperatures  since 
the  plate  is  mounted  in  an  inverted  position  near  the  top  of  the  test 
section  to  reduce  convection  currents  in  the  thermal  layer.  Occasional 
sanding  and  polishing  of  the  test  plate  was  done  to  remove  scale, 
a possible  source  of  hot  spots. 

D.  Typical  Run  Sequence 

For  a typical  run  the  following  steps  were  usually  taken:  With 

oscilloscopes  and  monitoring  instrumentation  adequately  warmed  up,  with 
test  plate  oven  temperature  in  the  desired  range,  the  plate  with  sting 
attached  is  removed  from  the  oven  and  (cold)  mounting  stud  affixed  to 
back  surface.  The  plate  assembly  is  then  brought  over  and  thrust  into 
the  test  section  and  fastened  to  its  ceiling  and  back  end,  while 
simultaneously,  final  adjustments  for  the  laser  Q-switching  are  made. 
The  driver  section  then  receives  its  final  charge  of  dry  Nj  for  the 

selected  pressure,  with  the  driven  section  at  atmospheric.  Simul- 
taneously, the  film  is  opened  in  the  darkened  test  section  building, 
and  recording  and  monitoring  systems  readied  for  the  run.  The  shot  is 
taken  usually  1-1/2  - 3 minutes  after  opening  of  the  oven. 
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E.  Data  Generation 


Records  from  a typical  run  are  shown  in  Figure  5.  The  shock 
passage  over  successive  gages  is  pictured  in  the  oscilloscope  traces 
and  permits  setting  of  the  delayed  triggering  of  the  laser  to  capture 
the  shock  at  the  test  section  window.  The  interferogram  for  this  run 
is  shown,  and  these  records  constitute  the  principal  raw  data  of  an 
experiment. 


A closer  look  at  the  interferogram  reveals  a number  of  interesting 
features  and  suggests  a picture  of  the  interaction  process.  The  un- 
disturbed region  and  the  uniform  region  behind  the  shock  in  the  cold 
air  are  characterized  by  the  (more  or  less)  straight  fringes.  The 
very  fine  and  closely  spaced  fringes  "above"  the  heated  plate  characterize 
the  heated  air  layer.  Last  but  not  least  is  an  interesting  feature  also 
observed  by  Griffith'1*  of  shock-thermal  layer  interaction:  one  notes  a 

"transition"  region,  a rounding,  shifting  of  fringes  through  the 
shocked  region  as  one  goes  from  the  base  of  the  shock  in  the  cold  region 
to  the  thermal  layer.  Some  of  cur  rounding  fringes  seem  broken  into 
straight  line  segments  and  suggest  that  the  primary  shock  degenerates 
into  a set  of  weak  shocks  or  pressure  waves  traversing  the  "transition" 
region  near  the  thermal  layer.  !i.  this  run  the  plate  (surface) 
temperature  was  determined  to  be  "05oK,  with  sound  speed  near  the 
plate  much  greater  than  shock  speed  in  the  cold  air;  thus  the  shock 
does  not  penetrate  to  the  plate  surface.  However,  a weak  disturbance 
is  to  be  noted  throughout  the  height  of  the  layer  if  one  follows 
individual  fringes  in  the  layer.  The  disturbance  apparently  is  a 
continuation  into  the  layer  of  the  leading  portion  of  the  foot  of  the 
shock  disturbance.  The  small  fringe  shifts  noted  in  the  layer 
represent  the  flow  property  adjustment  to  either  side  of  the  disturbance. 

A second  interferogram,  using  the  "infinite- fringe"  setting,  is  shown 
in  the  Figure  6.  A number  of  runs  were  made  using  this  setting  because 
of  greater  simplicity  in  reduction,  although  generally  less  accurate. 
However,  some  features  of  the  flow  are  exhibited  more  prominently 
than  in  the  regular  multifringe  setting  of  Figure  5.  Particularly  one 
notices  in  the  transition  zone  the  fanning  out  of  the  fringes  coming 
down  from  the  main  shock.  Since  in  the  infinite  fringe  setting  the 
fringes  represent  contours  of  constant  density,  one  has  more  explicitly 
the  picture  mentioned  above  of  the  degeneration  of  the  shock  into  a 
set  of  waves. 


F.  Data  Reduction 


The  (average)  shock  velocity  at  the  test  location  is  determined 
as  mentioned  previously  from  the  time  interval  between  shock  passage 
over  gage  stations  of  known  separation.  In  turn,  all  the  properties 
behind  the  shock  are  known  in  the  ideal  gas  case,  which  is  assumed 
throughout  the  present  work  for  the  weak  shocks  and  the  thermal  layer 
temperatures  involved. 
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Fringes  in  an  interferogram  are  "read  out"  on  a large  screen 
digitizing  comparator.  The  interferogram  reduction  requires  the 
fringe  shift  measurement  in  the  desired  region  of  flow  plus  pre- 
sumably known  initial  conditions  of  the  cold  air.  The  fringe  shift 
equation  relates  the  fringe  shift  to  the  density  change  through  the 
Gladstone-Dale  relation  and  may  be  written 


where 


S = 3 (L/A) 


6 = fringe  shift  from  undisturbed  fringe  (spacing) 

3 = (a  constant  for  the  gas  related  to  the 

refractivity  K)  = Kpg 

L = geometric  path  through  the  disturbance 
A a wavelength  of  source  light 


Pb'  pa'  ps  = Cities;  P^  is  density  producing  the 
fringe  shift  6 from  initial  density  p , 
with  p the  standard  density  at  273. 2°k  and 
S 5 

1.01325  x 10  Pa 


Use  of  the  ideal  gas  equation  of  state  p = pRT/M  leads  to  an  expression 
for  temperature  T in  terms  of  the  fringe  shift 


Ts/Tb  = ps/pb  tgL  + p T J * 

rs  a 


The  pressure  in  the  undisturbed  region  is  the  ambient;  that  in  the 
post-shock  region  is  determined  from  knowledge  of  the  shock  velocity; 
pressures  are  assumed  to  be  applied  down  through  the  thermal  layer  for 
these  uniform  regions.  For  the  multifringe  interferogram,  the  fringe 
shifts  from  undisturbed  fringe  in  both  regions  were  calculated  from  an 
average  slope  determined  for  the  undisturbed  fringe.  Ambient  tempera- 
ture before  disturbance  is  assumed  to  be  300°K, although  some  tempera- 
ture change  of  the  air  might  be  anticipated  in  the  test  section  from 
the  presence  of  the  heated  plate.  Ambient  temperature  in  the  post- 
shock region  is  again  from  shock  tables.  Thus  using  such  data  in  the 
temperature-fringe  shift  relation  allows  determination  of  the  tempera- 
ture profiles  in  pre-and  post-shock  regions  of  the  thermal  layer. 


The  region  at  the  base  of  the  chock  and  adjacent  to  the  thermal 
layer  requires  separate  consideration.  Using  the  interferogram  of 
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Figure  6 (with  the  infinite-fringe  setting)  for  discussion,  we  noted 
that  the  fringes  represent  lines  of  constant  density.  Thus  in  this 
situation,  an  inverse  analogy  is  suggested  to  the  Mach  waves  or 
characteristics  in  steady  two-dimensional  flow  over  a gradual  turn 
which  intersect  and  coalesce  into  the  main  shock  wave.  Here,  in 
the  flow  over  the  thermal  layer  - with  temperature/ sound  speed 
gradient  - we  see  dist  nuances  from  the  main  shock  reaching  down  into 
the  layer.  The  deep-lying  subsonic  portions  of  the  layer  must  affect 
the  disturbances  such  that  acoustic  waves  propagate  upstream  in  these 
portions,  as  evidenced  by  the  small  "breaks"  in  the  fringes.  Above 
these  portions  the  propagation  velocities  must  be  supersonic,  and  the 
observed  "fan"  of  weak  shock  waves  propagate.  Through  velocity 
continuity  considerations  across  any  layer,  one  expects  the  leading 
or  "precursor"  nature  of  the  shock  pattern.  Moreover,  a succeeding 
weak  wave  moves  at  slower  relative  velocity  (or  has  less  an  angle) 
since  the  wave  leading  has  heated  up,  compressed,  and  set  into  motion 
the  region  ahead  of  the  succeeding  wave;  hence  it  sees  a smaller 
driving  pressure  ratio.  In  sum,  the  interferogram  portrays  in  rather 
graphic  form  the  complex  shock-thermal  layei  interaction. 

In  view  of  this  picture  of  a "transition"  zone,  we  may  make  use 
of  the  weak  shock  relations  to  assign  temperatures  in  this  region.  An 
implicit  but  not  unreasonable  assumption  is  that  the  region  behind  each 
weak  shock  is  uniform.  In  the  deep  lying  layers  beneath  the  foot  of 
the  shock  the  pressure  applied  in  the  layers  must  vary  somehow  along 
the  plate  direction.  Hence,  adjustment  of  the  temperature  in  these 
regions,  from  the  uniform  pre-  and  post-shock  conditions,  is  unresolved 
in  this  work.  The  density  "jump"  at  a fixed  height  in  the  deep-lying 
layers  could  be  gotten  from  the  fringe  shift,  as  usual;  but  this  has 
not  been  done.  For  the  region  spanned  by  the  fan  of  weak  shocks,  the 
weak  shock  relations  are 


Ap_  _ 1_  A£ 
P„  = Y P„ 


CY-l) 


Ap 


Y P 

a ra 

which  coupled  with  the  fringe  shift  relation 


Ap 

P 


s X 


P 3L 
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give  the  working  equations 


[(Y-l) 


— "6  + 1]  T 
p BL  1 a 
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and 


A . 

= Pc  oT  5 + p„  • 

D S pL  Si 

The  first  calculation  begins  from  the  known  undisturbed  conditions. 
Succeeding  flow  values  are  obtained  by  iteration.  The  multifringe 
interferograms  are  treated  in  a similar  manner.  The  results  of  such 
manipulations  and  of  the  experiments  follow. 


III.  RESULTS  AND  DISCUSSION 

Shown  in  Figure  7 are  examples  of  completely  reduced  data  runs. 
Accompanying  each  is  the  interferogram  for  the  run.  The  shock  velocity 
and  Mach  number  are  noted  in  the  figures.  The  fringes  sketched  are 
scaled  4:1  of  actual  physical  dimensions,  as  determined  from  the 
photographs  and  the  known  heated  plate  thickness,  better  to  exhibit 
the  temperature  data.  Temperature  profiles  in  pre-  and  post-shock 
regions  are  illustrated.  The  temperature  dependence  for  both  regions 
appears  generally  to  be  an  exponential  decrease  with  increasing  height 
above  the  plate,  for  about  the  first  0.75  cm,  as  the  plots  of  Figure  8 
show,  and  thereafter  a decrease  at  a much  slower  rate  until  the  ambient 
temperature  is  reached.  Of  note  is  the  relatively  good  agreement  of 
heated  plate  surface  temperature  T^  to  either  side  of  the  shock.  The 

infinite  fringe  interferograms  exhibited  greater  disparity  in  the  T^, 

but  values  were  well  within  3%  of  each  other.  For  the  multifringe 
interferogram  of  Figure  7B,  one  may  draw  in  local  temperature  contours 
in  the  transition  zone,  and  most  tend  to  straight  lines,  corroborating 
features  of  the  infinite  fringe  interferograms.  In  the  transition 
zone,  in  a majority  of  the  runs  the  final  temperature  reached  dis- 
agrees with  the  post-shock  region  temperature  predicted  from  shock 
tables  based  on  the  measured  shock  velocity.  Some  of  the  discrepancy 
may  be  attributed  to  the  cumulative  error  in  stepping  through  the 
region.  Other  possible  sources  for  error  are  discussed  below. 

Accuracy  of  the  numbers  generated  is  limited  by  several  factors. 
Our  knowledge  of  the  initial,  undisturbed,  ambient  conditions, 
particularly  with  heated  plate  in  position,  is  a question.  Beyond 
this,  the  shock  speed  is  accurate  to  ± 1%  as  mentioned  earlier; 
and  comparator  readings  for  fringe  coordinates  appear  to  be  accurate 
and  reproducible.  The  readings  appeared  to  represent  the  fringe 
shapes  faithfully,  so  that  extrapolation  of  undisturbed  (straight) 
fringes  for  determination  of  fringe  shift  should  suffer  no  more  than 
is  usual.  Of  greater  impact  here  is  actually  some  curvature  of 
fringes  in  the  supposedly  uniform  regions,  a feature  noted  in  previous 
work  in  this  facility1*  but  not  completely  explained  away.  Areas  with 


11.  J.  H.  Spurk  and  J.  M.  Bartos , "Interferometric  Measurement  of  the 
Nonequilibrium  Flow  Over  a Cone."  Phys.  Fluids  £,  1278-85  (1966). 
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noticeable  fringe  curvature  were  generally  avoided,  particularly  at 
the  field  extremes. 

Significant  for  the  infinite  fringe  interferograms  was  that  the 
interferometer  setting,  as  it  turned  out  after  the  run  on  several 
occasions,  had  drifted  very  slightly,  with  the  result  that  several 
fringes,  or  some  peculiar,  large  fringe  pattern  had  set  in,  rather 
than  the  more  or  less  evenly  illuminated  field  desired  (see  Figure  6 
and  Figure  7C,  post-shock  region).  Thus,  the  "zero"  reference  fringe 
from  which  the-  ensuing  visible  fringes  have  shifted  is  a somewhat 
questionable  quantity.  In  order  to  reduce  the  data  from  such  runs, 
a bit  of  interpretation  was  used  to  guess  the  reference  fringe,  based 
on  the  expected  fringe  shift  through  the  shock,  the  height  of  the 
thermal  layer  above  the  heated  plate,  plus  a count  of  fringes  through 
the  transition  zone.  Notwithstanding  these  uncertainties,  one  notes 
disturbances  behind  the  shock  in  all  the  interferograms.  One  dis- 
turbance was  traced  to  a shock  reflection  from  a protruding  gasket, 
subsequently  re-positioned,  upstream  at  the  test  section  joint.  A 
persisting  disturbance  was  the  shock  reflection  off  the  thickened 
"leading  edge"  of  the  plate  caused  by  the  thermal  layer  air  currents 
around  the  edge.  This  effect  is  illustrated  in  the  Figure  9. 

(Observation  of  the  shocked  flow  over  the  cold  plate  itself  showed  no 
(blunt)  leading  edge  effect  at  the  test  station,  centered  25cm  down- 
stream of  the  leading  edge.)  Clearly,  for  very  accurate  work  attention 
to  the  above  mentioned  uncertainties  plus  others  would  be  required 
either  to  minimize  or  to  estimate  their  effect  on  the  experiment.  The 
present  work  lias,  however,  been  exploratory  and  demonstrates  some  of 
the  possibilities  and  the  pitfalls  associated  with  the  shock/thermal 
layer  interaction  experiment  and  the  interferometric  technique  for  the 
desired  temperature  measurement.  In  the  main,  because  of  its  familiarity 
and  relative  simplicity  in  data  reduction  procedures,  the  technique 
seems  well  suited  for  these  studies. 

Finally,  some  observations  of  the  connection  of  the  present 
experiments  and  results  with  field  results  should  be  made.  The 
results  and  features  protrayed  in  this  work  bear  just  some  resemblance 
to  precursor  behavior  in  the  field  as  depicted  in  Figure  1.  Also,  we 
have  not  attempted  a comparison  with  the  relation  between  precursor 
angle  and  temperature  ratio,  because  of  our  non-uniform  thermal  layer. 
However,  some  tenuous  correspondences  might  be  drawn:  The  pressure 

levels  in  this  work  are  appropriate  to  nuclear  blasts  at  GR  'v  600  m 
and  beyond,  where  precursor  behavior  seems  less  pronounced.  In  Figure  10 
are  examples  of  pressure  traces,  from  a shot  with  precursor  formation, 
for  such  distances,*  which  are  described  as  typical  "clean-up"  wave- 
forms as  the  shock  moves  farther  and  farther  from  the  source. 

The  ground  level  records  OB  and  2B  show  a rather  gradual  rise  to 
a plateau,  followed  by  some  indication  of  a cusp  (more  pronounced  in 
2B)  (and  subsequent  drop-off  associated  with  the  spherical  blast). 

The  gradual  pressure  rise  might  correspond  to  the  undetermined 
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pressure  development  implied  by  the  interf erograms  in  passing  from  the 
pre-shock  region  to  the  post-shock  region  along  the  heated  plate  sur- 
face. The  cusp  might  then  be  identified  with  the  small  break  mentioned 
earlier  in  the  fringes  deep  in  the  heated  layer. 

The  above  ground  stations  2B10  and  2B50  show  the  pressure  rise 
and  increase  in  pressure  to  a second  sharp  rise  before  the  pressure 
decrease  of  the  spherical  blast.  The  transition  zone  behavior  depicted 
in  the  interf erograms  would  correspond  thusly:  the  initial  and  rising 

pressures  are  those  associated  with  the  first  and  subsequent  waves  of 
the  transition  region.  The  sharper  second  pressure  rise  at  greater 
height  would  signify  the  lessened  influence  at  greater  height  of  the 
thermal  laye:  i.e.,the  lessened  extent  of  the  wave  "fan"  from  the 

main  shock  and  the  more  normal-shock  behavior  of  the  disturbance. 

Additionally,  shock  contours  plotted  from  direct  shock  photography 
of  nuclear  bursts  show  some  interesting  features.  Figure  11  reproduces 
a few  pertinent  ones  from  a Figure  in  Reference  1,  which  shows,  for  a 
shot,  shock  contours  and  their  heights  at  various  ground  ranges;  times 
after  burst  are  also  indicated.  In  particular,  we  note  the  multiple 
shocks  fanning  down  from  single,  main  shocks,  bearing  rather  clear 
similarity  to  the  shocks  of  the  transition  zone  in  the  present  experi- 
ments. In  the  field  tests  a number  of  density  gradients  were  seen 
behind  the  leading  precursor  wavefront  at  any  instant.  Many  of  these 
were  thought  to  be  local  disturbances.  The  results  of  our  experiments 
suggest  another  likelihood  also:  that  the  density  gradients,  or 
equivalently,  the  thin  regions  constituting  the  transition  zone  are 
intimately  connected  with  the  interactions  between  the  main  shock  and 
the  thermal  layer  and  not  solely  with  particular,  local  disturbances. 

We  stress  that  these  remarks  are  somewhat  speculative,  and  no  additional 
evidence  is  presently  in  hand  to  support  them. 


IV.  CONCLUSIONS 

This  report,  based  on  observations  using  a Mach  Zehnder  inter- 
ferometer, examines  some  features  of  the  interaction  of  a plane  shock 
generated  by  a shock  tube  with  a thermal  layer  created  over  a heated 
plate.  The  technique  is  similar  to  that  of  Griffith  who  worked  at 
much  lower  temperatures.  Temperature  profiles  through  the  thermal 
layer  are  determined  in  the  pre-  and  post-shock  regions.  The  inter- 
ferograms  show  a feature  of  the  interaction  - a "transition  zone"  which 
is  interpreted  as  a decomposition  of  the  main  shock  into  a fan  of 
weak  shock  or  Mach  waves  reaching  down  into  the  thermal  layer.  Some 
of  the  waves  "lead"  the  main  shock  and  would  be  associated  with  the 
normally  observed  "precursor"  shock.  Some  assumptions  are  made  to 
assign  temperatures  in  this  transition  zone.  Tentative  connections 
of  the  experimental  results  with  field  results  are  discussed. 
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